#include "Eigen/Dense"
#include <iostream>

int main() {
    // 示例1：基础矩阵运算
    Eigen::Matrix2d mat;
    mat << 1, 2,
           3, 4;
    
    Eigen::Vector2d vec(5, 6);
    Eigen::Vector2d result = mat * vec;

    std::cout << "矩阵乘法结果：\n" << result << std::endl;

    // 示例2：求解线性方程组
    Eigen::Matrix3f A;
    Eigen::Vector3f b;
    A << 1, 2, 3,
         4, 5, 6,
         7, 8, 10;
    b << 3, 3, 4;

    Eigen::Vector3f x = A.fullPivLu().solve(b);

    std::cout << "线性方程组的解：\n" << x << std::endl;  

    return 0;
}